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• Premise of the study: Rhizopogon kretzerae and R. salebrosus (Rhizopogonaceae) are ectomycorrhizal fungi symbiotic with 
pines and the mycoheterotrophic plant Pterospora andromedea (Ericaceae). Microsatellite loci will allow population genetic 
study of fungal hosts to P. andromedea. 

• Methods and Results: Shotgun pyrosequencing of R. kretzerae DNA resulted in primer development of 23 perfect microsatel- 
lite loci and screened across two populations each for R. kretzerae and R. salebrosus. Twelve loci were polymorphic in 
R. kretzerae populations, and 1 1 loci cross-amplified in R. salebrosus populations. For R. kretzerae and R. salebrosus, number 
of alleles was one to eight and one to nine, respectively, and observed heterozygosity ranged from 0.00-0.57 and 0.00-0.70, 
respectively. 

• Conclusions: These are the first microsatellite loci developed for any species within Rhizopogon subgenus Amylopogon. These 
microsatellite loci will be used in conservation genetic studies of rare to endangered eastern populations and to compare plant 
and fungal population genetic structure at different hierarchical levels. 

Key words: conservation genetics; GS-FLX Titanium; micro satellites; rare fungi; Rhizopogon kretzerae; Rhizopogon 
subgenus Amylopogon. 



Rhizopogon Fr. (Rhizopogonaceae) is a fungal genus of spe- 
cies that form ectomycorrhizal symbioses with conifers and 
are pivotal components of temperate forest ecosystems (Molina 
et al., 1999). Four species (R. kretzerae Grubisha, Dowie & 
Mill., R. salebrosus A. H. Sm., R. ellenae A. H. Sm., and R. 
arctostaphyli A. H. Sm.) in Rhizopogon subgenus Amylopogon 
(A. H. Sm.) Grubisha & Trappe are obligately ectomycorrhizal 
with several pine species and are the symbiotic mycobiont host 
to the mycoheterotrophic plant, Pterospora andromedea Nutt. 
(Ericaceae; Bidartondo and Bruns, 2002; Dowie et al., 2011, 
2012; Hazard et al, 2012; Grubisha et al, 2014), a North Amer- 
ican endemic divided into eastern and western populations 
(Bakshi, 1959). Pterospora andromedea primarily associates 
with two closely related Rhizopogon species: western P. andro- 
medea are symbiotic primarily with R. salebrosus while eastern 
populations are found with R. kretzerae (Grubisha et al., 2014). 
Rhizopogon kretzerae is currently only known from collections 
of P. andromedea roots and only found in association with 
Pinus strobus L. (Hazard et al., 2012; Grubisha et al., 2014). 
Eastern P. andromedea populations have undergone recent de- 
creases in population size and possible local extinction from a 
variety of anthropogenic factors (Schori, 2002). 

1 Manuscript received 22 March 2014; revision accepted 20 May 2014. 
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Foundation (DEB- 10503 15 to M.R.K. and DEB- 1050292 to S.L.M.). 
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Copopulation genetic studies of the fungal host and myco- 
heterotrophic plant will (1) provide information about levels of 
genetic variation and population structure at different hierarchi- 
cal levels in plant and fungal populations, and (2) be important 
in developing conservation management plans for the rare to 
endangered eastern populations of both species. While micro- 
satellite loci have been previously developed for distantly re- 
lated species in three of the five Rhizopogon subgenera (Kretzer 
et al., 2000, 2004; Grubisha et al, 2005), the microsatellite loci 
reported here, which were developed from R. kretzerae and R. 
salebrosus, are the first for any species in Rhizopogon subgenus 
Amylopogon. 

METHODS AND RESULTS 

Genomic DNA (gDNA) was isolated from 60 mg of silica gel-dried pure 
cultures of R. kretzerae isolate Miller207 (Michigan; Table 1) and R. salebro- 
sus isolate Miller322 (Washington; 46.20545°N, 117.7704°W) separately 
using the DNeasy Plant Mini Kit (QIAGEN, Valencia, California, USA). These 
two cultures were obtained from P. andromedea monotropoid roots following 
Grubisha et al. (2014). Each culture was ground in liquid nitrogen to a fine 
powder. Twice the recommended volume of buffer API, RNaseA, and buffer 
AP2 were used. Library construction and 454 shotgun pyrosequencing were 
performed by the Advanced Genetic Technologies Centre (AGTC) at the Uni- 
versity of Kentucky (Lexington, Kentucky, USA). Libraries for each species 
were prepared using the GS FLX+/XL+ Titanium Rapid Library Preparation Kit 
(454 Life Sciences, a Roche Company, Branford, Connecticut, USA). Instead 
of nebulizing, gDNA was sheared using a Bioruptor NGS (UCD-600TS; 
Diagenode, Denville, New Jersey, USA) by sonicating at 4°C with six cycles of 
5 s on and 90 s off. Shotgun 454 pyrosequencing was performed on 25% of a 
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plate on a Titanium GS-FLX (454 Life Sciences, a Roche Company) yielding a 
total of 152,007 reads with an average length of 414 bp for R. kretzerae and a 
total of 92,252 reads with an average length of 517 bp for R. salebrosus. At the 
AGTC, raw sequence reads were initially processed using prinseq-lite.pl (http:// 
prinseq.sourceforge.net/) and involved removing all (1) reads with >2% length 
with uncalled bases (N), (2) duplicate reads, (3) reads with Q < 17, (4) reads too 
short in length [average length - 2(standard deviation)], and (5) reads too long 
in length [average length + 2(standard deviation)]. Processed reads were as- 
sembled de novo using Newbler version 2.8 (454 Life Sciences, a Roche Com- 
pany) at the AGTC, which produced 10,190 contigs for R. kretzerae (average 
contig size = 1502 bp, N50 = 1701 bp) and 5178 contigs for/?, salebrosus (aver- 
age contig size = 1302 bp, N50 = 1428 bp). 

Perfect microsatellite repeats were identified from the assembled contigs using 
MS ATCOMMANDER version 1 .0.8 (Faircloth, 2008) for R. kretzerae and R. sale- 
brosus, yielding 43 and 21 dinucleotides with at least six repeat units, 137 and 48 
trinucleotides with at least five repeat units, 34 and 12 tetranucleotides with at least 
four repeat units, five and 1 1 pentanucleotides with at least four repeat units, and 
five and 1 1 hexanucleotides with at least four repeat units, respectively. Within 
MS ATCOMMANDER, Primer3 (Rozen and Skaletsky, 2000) was used to gener- 
ate primers for all loci using default parameters except that difference in acceptable 
primer melting temperature was 2°C with an optimum of 60°C; maximum poly X 
was set to three; GC clamp = yes; and a target PCR product size of 90-210 bp. Se- 
lection of loci for initial screening was based on Primer3 output tabulated in a 
spreadsheet by MSATCOMMANDER, and those that had no lowercase letters in 
the primer sequences that would reflect nucleotide mismatches across reads in as- 
sembled contigs. Forty-eight loci (38 identified from R. kretzerae and 10 from R. 
salebrosus) were initially screened across four samples (two R. kretzerae, two R. 
salebrosus; Table 1) for positive PCR products using agarose gel electrophoresis. 
Loci were PCR amplified using 0.1 jiM of each forward and reverse primer in a 
10-jiL reaction volume that included lx CoralLoad PCR Buffer with 1.5 mM 
MgCl 2 (QIAGEN), 200 \iM each dNTP, 2.5 units Taq DNA Polymerase 
(QIAGEN), and 1.0 jiL of 1 : 10 diluted genomic DNA. Thermocycler parameters 
were 94°C for 3 min, 35 cycles of 94°C for 30 s, 55°C for 30 s, and 72°C for 30 s, 
followed by a final extension for 10 min at 72°C. PCR products were visualized in 
a 2% agarose gel stained with GelRed (Phenix Research Products, Candler, North 
Carolina, USA) on a UV transilluminator. Of the 48 loci assessed, 27 loci (23 from 
R. kretzerae and four from R. salebrosus) yielded no more than one or two bright 
PCR product bands within the approximate expected size range (target size + 40 bp, 
-10 bp) and no PCR products outside the expected size range. The lower threshold 
was set to target size -10 bp due to the relatively low number of repeats. 

These 27 loci were screened in 65 isolates from two populations each for R. 
kretzerae (n = 14 and n = 11; Table 1) and R. salebrosus (n = 20 both; Table 1), 
and for cross-amplification in two distantly related Rhizopogon subgenus Amylo- 
pogon species that are also hosts for P. andromedea: R. arctostaphyli (n = 2 for 
both; Table 1) and R. ellenae (n = 2 for both; Table 1). Pterospora andromedea 
ectomycorrhizal roots were collected and fungal DNA was isolated as previously 
described (Dowie et al., 2012; Grubisha et al., 2014). Due to the extremely small 
sample of rootball collected, voucher specimens in herbaria were not made except 
for samples from Michigan (DAOM 242738) and Quebec (DAOM 242739) that 
were previously deposited at the Canadian National Mycological Herbarium 
(Grubisha et al., 2014). The forward primer was 5' end-labeled with one of four 
dyes in the ABI DS-33 dye set (Table 2). Loci were PCR amplified in multiplex 
reactions using the QIAGEN Multiplex PCR Kit in a 10-jiL volume with lx 
QIAGEN Multiplex PCR Master Mix, 50 nM each primer, 0.25-0.5x Q-solution, 
and 1.0 JiL of 1 : 10 diluted genomic DNA. Touchdown thermocycler conditions 
were: 95°C for 15 min; 10 cycles of 94°C for 30 s, 67°C for 90 s, decreasing 1°C 
each cycle, and 72°C for 30 s; 25 cycles of 94°C for 30 s, 57°C for 90 s, 72°C for 
30 s; with a final extension of 45 min at 60°C. Primers with a higher or lower 
touchdown annealing profile are listed in Table 2. Fragment analysis was con- 
ducted using the GeneScan 500 LIZ Size Standard (Applied Biosystems, Foster 
City, California, USA) on an ABI 3730 DNA Analyzer (Applied Biosystems) by 
the Biotechnology Resource Center (BRC) at Cornell University. Allele sizes 
were called manually using the Microsatellite Plugin in Geneious version R6. 1 .5 
(Drummond et al., 201 1). Loci selected for the final set were PCR amplified inde- 
pendently to verify allele size ranges. 

Of the 23 loci developed in R. kretzerae and screened with fluorescently la- 
beled primers in R. kretzerae populations, 12 loci were polymorphic (Table 1), 
seven loci were monomorphic, and four loci had excessive stutter. In R. sale- 
brosus populations, 21 of the 23 loci developed in R. kretzerae successfully 
PCR amplified and 1 1 loci were polymorphic, four of which were monomor- 
phic in R. kretzerae populations (Table 1). The four loci isolated in R. salebro- 
sus did not PCR amplify in R. kretzerae populations screened. In R. salebrosus 
populations, only one locus (Rsal2) was polymorphic (Table 1). A few loci 
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Table 2. Characteristics of 12 polymorphic and five monomorphic micro satellite loci developed in Rhizopogon kretzerae and one polymorphic locus 
developed in R. salebrosus. 

GenBank 



Locus a 


5' end-labeled dye 




Primer sequences (5'-3') 


Repeat motif 


Allele size range (bp) 


T a (°Q 


accession no. 


Rkre3 


NED 


F 
R 


TACCTTGCCCACTATCGACG 
TAATCGATGGACCCGCTGAG 


(TGC) 5 


200-221 


60 


KJ443650 


Rkre5 


VIC 


F 

R 


AGACCAATACGCGTGAACAC 
CAGAGCAATCCGAACACGAG 


(CG) 6 


179-183 


57 


KJ443651 


Rkre6 


PET 


F 

R 


CCGTGACAATCCAACTGCAC 
GGCAACACTGAATCTGGAGC 


(ACG) 6 


218-221 


57 


KJ443652 


Rkrel4 


VIC 


F 

R 


AATGCTGGACTTCTTGCACG 
TGGGCATCATGGTGTAGGAC 


(AG) 6 


192-196 


60 


KJ443653 


Rkrel5 


6-FAM 


F 

R 


TGTTTGGTCGTGATCTTGTAGC 
ACCTTCGTAGTCCAGTGCTG 


(AT) 7 


169-181 


52 


KJ443654 


Rkrel7 


PET 


F 
R 


AAGCAATAAAGGAATCTGCAGC 
CCGGATGTCTAACACTGGTC 


CAT) n 


176 


57 


KJ443659 


Rkrel8 


6-FAM 


F 
R 


TCAATGACCGGTCCACTCTC 
GTCAGTAGGTGTATTCTTGCTGG 


(ATGC) 5 


120-128 


60 


KJ443655 


Rkrel9 


PET 


F 
R 


GATGCAAGGCTCGTTGTCTC 
GCAAGTTCTGGAGCAAGGAC 


(TCT) 10 


202 


57 


KJ443660 


Rkre20 


VIC 


F 
R 


TGCTTCT C AAAC AAC C G AGG 
GGTGGTAGCGATTGAGGGAC 


(ACC) 6 


114 


57 


KJ443661 


Rkre23 


PET 


F 
R 


ACTTTGATGCAGGCGAGAAC 
TGACCGTCACAAATCGTTCC 


(AGC) 5 


129-132 


57 


KJ443656 


Rkre24 


NED 


F 
R 


TCAGGGTTGGACTCGTGAAG 
CGCGCCTGAAGTTTGGAC 


(AGC) 7 


117-132 


52 


KJ443657 


Rkre27 


6-FAM 


F 
R 


ATGGAGGATGAGAGTGCGG 
CGTCGTGACATGAGTGCTTG 


(ACG) 6 


108-111 


57 


KJ443658 


Rkre29 


PET 


F 
R 


AGGTTGGCTGTGCATGAATG 
ACACTAGCCTGCTTCCTCAG 


(AT) 6 


188 


57 


KJ443662 


Rkre34 


PET 


F 
R 


TATAGACTCGGACGCACACC 
AGGTTCCCGTGAGTGATATCG 


(CA) 10 


119-131 


57 


KJ443663 


Rkre36 


NED 


F 
R 


AAGCAGGGTACAGACTCACC 
TCTGGATCGTTTGCAATGCC 


<TC)„ 


173-215 


57 


KJ443664 


Rkre37 


VIC 


F 

R 


TGGTGGTGGGTATGAGTTGG 
TTTGTCAAAGCAGGCTACGC 


(TGG) 9 


183-217 


52 


KJ443665 


Rkre38 


NED 


F 

R 


ATGCAGGTAGACGGGACTTG 
TACTCGTCGCTACAGTCCG 


(GTG) 6 


209 


57 


KJ443666 


Rsal2 


NED 


F 

R 


TGTCCTCAAGAGATCCGCAG 
AGATTCATCGAAACCGCACG 


(TCG) 5 


110-119 


57 


KJ443667 



Note: T a = annealing temperature. 

a Rkre = isolated from R. kretzerae; Rsal = isolated from R. salebrosus. 



cross-amplified in R. arctostaphyli and R. ellenae isolates and may prove to be 
useful when screened at a larger scale (Table 1). Locus Rkrel9 had stutter prob- 
lems in both R. kretzerae and R. salebrosus but was polymorphic in R. arcto- 
staphyli and R. ellenae (Table 1). 

Number and size range of alleles, and observed (H 0 ) and expected heterozygosi- 
ties (H e ) were calculated in GenAlEx version 6.5 (Peakall and Smouse, 2006, 2012; 
Table 1). When both populations for each species were considered (R. kretzerae, 
n = 25; R. salebrosus, n = 40), for the loci isolated in R. kretzerae (Rkre3-Rkre38; 
Table 1), there were one to eight alleles, and H e ranged from 0.00 to 0.64 in 
R. kretzerae populations (Table 1); and there were one to nine alleles and H e values 
of 0.00 to 0.76 in R. salebrosus populations (Table 1). For locus Rsal2 there were 
one to three alleles, and H e ranged from 0.00 to 0.10. Scoring errors due to stutter- 
ing, large allele dropout, and null alleles were checked using MICROCHECKER 
2.2.3 (van Oosterhout et al., 2004). For R. kretzerae populations, two loci showed 
signs of possible null alleles (Rkre5 in Quebec and Rkre37 in Michigan). In the 
combined R. salebrosus populations, possible null alleles and scoring errors due to 
stuttering were detected at four loci (Rkre3, Rkrel7, Rkre27, Rkre34). 



CONCLUSIONS 

These are the first microsatellite loci developed for any spe- 
cies in Rhizopogon subgenus Amylopogon. Presently we are us- 
ing these loci to assess the conservation genetic status of R. 
kretzerae, which will be important for developing conservation 



management plans for eastern R andromedea. We are also con- 
ducting copopulation genetic studies to determine if this fungal 
host and mycoheterotrophic plant share similar patterns of pop- 
ulation genetic structure and variation at different hierarchical 
levels. Cross-amplification of loci isolated in R. kretzerae with 
R. salebrosus and two distantly related Rhizopogon subgenus 
Amylopogon species suggests they may be useful in studies of 
other species in this subgenus. 
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